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Abstract. File-sharing networks are distributed systems used to disseminate 
files among nodes of a communication network. The general simple principle 
of these systems is that once a node has retrieved a file, it may become a server 
for this file. In this paper, the capacity of these networks is analyzed with a 
stochastic model when there is a constant flow of incoming requests for a given 
file. It is shown that the problem can be solved by analyzing the asymptotic 
behavior of a class of interacting branching processes. Several results of in- 
dependent interest concerning these branching processes are derived and then 
used to study the file-sharing systems. 



File-sharing networks are distributed systems used to disseminate information 
among some nodes of a communication network. The general simple principle is 
the following: once a node has retrieved a file it becomes a server for this file. An 
improved version of this principle consists in splitting the original file into several 
pieces (called "chunks") so that a given node can retrieve simultaneously several 
chunks of the same file from different servers. In this case, the rate to get a given 
file may thus increase significantly as well as the global capacity of the file-sharing 
system: a node becomes a server of a chunk as soon as it has retrieved it and not 
only when it has the whole file. These schemes disseminate information efficiently 
as long as the number of servers is growing rapidly. 

The paper investigates the maximal throughput of these file-sharing networks, 
i.e., if the system can cope or not with a constant flow of incoming requests for 
a given file. Two cases are considered, either the file consists of one chunk or the 
file is split in two chunks which are retrieved sequentially, which we will refer to 
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as a linear file-sharing network. It is assumed that arrival times are Poisson and 
chunk transmission times are exponentially distributed. In this setting, the stability 
property of the file-sharing system is expressed as the ergodicity property of the 
associated Markov process. Even in this simple framework, mathematical studies 
are quite scarce, see Qiu and Srikant [14] , Simatos et al. [18], Susitaival et al. [19] 
and references therein. 

The main technical difficulty in proving stability /instability results for this class 
of stochastic networks is that, except for the input, the Markov process has un- 
bounded jump rates, in fact proportional to one of the coordinates of the current 
state. If x servers have a chunk, since each of them can deliver this chunk, they 
globally provide this chunk at a rate proportional to x. For this reason, the classical 
tools related to fluid limits cannot be used easily in this setting. See Bramson [4], 
Chen and Yao [5] and Robert [15] for example. However, this class of processes 
is close to another important class of stochastic processes, namely branching pro- 
cesses, where a population of size x evolves at rate proportional to x. As it will 
be seen, a file sharing system distributing a file split in several chunks can be rep- 
resented as a Markov process associated to multi-type branching processes with 
interaction. 

Interacting Branching Processes. Consider the case of a file split into two 
chunks, one considers the case when a new user arriving in the network requests 
first chunk number 1 and then chunk number 2. We further assume that a server of 
type-1, i.e., having chunk 1 only, distributes it at rate [1%, while a server of type-2, 
having chunks 1 and 2, distributes only chunk 2 at rate ^2- For i — 0, 1, 2, let 
Xi(t) be the number of type-? servers at time t. A type-0 server is simply a user 
without any chunk. The crucial observation is that as long as there are requests 
for chunk i G {1,2}, then (Xi(t)) grows similarly as a branching process where 
individuals give birth to one child at rate fii, since each server offers a capacity fXi 
for the requested chunk. On the other hand, a new arrival in Xi corresponds to 
a departure in X\, since the new type-2 server was a type-1 server, thus deaths in 
X\ are governed by births in Xi. 



HlXi 
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Figure 1. Transition rates outside the boundaries of the file- 
sharing network with two chunks. 

The file-sharing network under consideration can thus be seen as a system of 
interacting branching processes where the births and the deaths of individuals are 
correlated. In the simpler case of one chunk, a description as a branching process 
(without interaction) has been used to analyze the transient behavior of these file- 
sharing systems. See Yang and de Veciana [3T], Dang et al. [B] and Simatos et 
al. [TS]. 

Sections [5] and |3] present results of independent interest concerning branching 
processes where individuals are killed at the instants (a n ) of a random point pro- 
cess. In Section [2] a criterion for the extinction of the branching process is obtained 
in terms of the sequence (a n ) and an asymptotic result is derived in this case. Sec- 
tion [3] considers the case where (cr„) is the sequence of birth instants of another 
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independent branching process. Several useful estimates are derived in this setting. 
These results are used to establish the stability results concerning file-sharing sys- 
tems with two chunks. The stability properties of a network with a single-chunk 
file are analyzed in detail in Section |4j The case of file-sharing networks with two 
chunks is detailed in Section [5] 

Acknowledgements. 

This paper has benefited from various interesting discussions with S. Borst, I. Nor- 
ros, R. Niinez-Queija, B.J. Prabhu, and H. Rcittu. 

2. Yule Processes with Killing 

A Yule process (Y(t)) with rate /x > is a Markovian branching process with 
Q-matrix 

(1) Qy{x, x + 1) = fix, \fx > 0. 

A Yule process is simply a pure birth process, where each individual gives birth to 
a child at rate \x. 

Yule processes are the basic branching processes that appear in analyzing the 
two-chunk network of Section [5] Actually a variant of this stochastic model will 
be needed, where some individuals are killed. In this section we study this model 
when killings are given by an exogenous process and occur at fixed (random or 
deterministic) epochs; in Section [3] killings result from the interaction with another 
branching process. 

In terms of branching processes, this killing procedure amounts to prune the tree, 
i.e., to cut some edges of the tree, and the subtree attached to it. This procedure 
is fairly common for branching processes, in the Crump-Mode- Jagers model for 
example, see Kingman [7]. See also Neveu [T2] or Aldous and Pitman pQ. 

A Yule Process Killed at Fixed Instants. Until the end of this section, it is 
assumed that, provided that it is non-empty, at epochs a n , n > 1, an individual 
is removed from the population of an ordinary Yule process (Y(t)) with rate nw 
starting with Y(0) = w G N individuals. It is assumed that (a n ) is some fixed 
non-decreasing sequence. It will be shown that the process (Wit)) obtained by 
killing one individual of (Y{t)) at each of the successive instants (er„) survives with 
positive probability when the series with general term (exp(— Hw^n)) converges. 
We denote 



where, for 1 < i < k and t > o",, Xj(t) is the total number of children at time t in 
the original Yule process of the ith individual killed at time cr^. In terms of trees, 
(W(t)) can be seen as a subtree of (Y(t)): for 1 < i < k, (Xi(t)) is the subtree of 
(Y(t)) associated with the zth particle killed at time cr;. 

It is easily checked that (Xi(t — cn),t > ai) is a Yule process starting with one 
individual and, since a killed individual cannot have one of his descendants killed, 
that the processes 



K = inf{n > 1 : W(a n ) = 0}. 
The process (W(t)) can be represented by 



n 



(2) 




{Xi(t)) = (Xi(t + (Ti),t>0), l<i<K, 
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are independent Yule processes. 

For any process (U(t)), one denotes 

(3) (M v (t)) d = (e-^*U(t)) . 

If (X(t)) is a Yule process with rate /iw, the martingale (M^(t)) converges almost 
surely and in L2 to a random variable M^{oo) with an exponential distribution 
with mean X(0), and by Doob's Inequality 

E(supM^(t) 2 ) < 2supE(M jf (f) 2 ) < +00. 

\t>0 J t>0 

See Athreya and Ney [5J- Consequently 

K 

e-» wt W{t) = My it) - e-wiM^t - ^)l {(7 .< t} , 
i=i 

and for any t > 0, 

re re 

l= i ' <= i ' 

Assume now that J2i>i e~^ w(Ti < +00: then the last expression is integrable, and 
Lebesgue's Theorem implies that (My/{t)) = (exp(— (iwt)W(t)) converges almost 
surely and in Li to 

K 

M w (oo) = M Y (00) -^Te-w^Mx.ioo). 
»=i 

Clearly, for some w* large enough and then for any w > w* , one has 

+00 

E w (M w (oo)) > w - e _w<r< > 0, 

i=l 

in particular P w (Myp(oo) > 0) > and P U) (W(t) > l,Vt > 0) > 0. If Y(0) = iu < 
w* and <ti > 0, then V w (Y(ai) > w* + 1) > and therefore, by translation at time 
(Ti, the same conclusion holds when the sequence (exp(— jAwCi)) has a finite sum. 
The following proposition has thus been proved. 

Proposition 2.1. Let (W(t)) be a process growing as a Yule process with rate /iw 
and for which individuals are killed at non- decreasing instants (o~ n ) with o\ > 0. If 

+00 

^e-""'"' < +00, 
i=i 

then as t gets large, and for any w > 1, the variable (exp(— iiyyt)W(t)) con- 
verges ¥ w -almost surely and in L2 to a finite random variable Mw(co) such that 
¥ w (M w (oo) > 0) > 0. 

The previous proposition establishes the minimal results needed in Section [5] 
However, Kolmogorov's Three-Series, see Williams |20j . can be used in conjunction 
with Fatou's Lemma to show that (W(t)) dies out almost surely when the series 
with general term (exp(— nwo~ n )) diverges. 
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3. Interacting Branching Processes 

In this section we study the Yule process of the previous section when killing 
times correspond to the birth times of some other branching process. The other 
branching process can be seen as a renewing Bellman-Harris process. 

Renewing Bellman-Harris process. In the rest of this section, nz,v > are 
fixed and {Z(t)) is a birth-and-death process whose Q-matrix Qz is given by 

(4) qz{z 1 z + 1) = [iz{z V 1) and qz(z, z — 1) = vz. 

In the rest of the paper n V m denotes max(rt, m) for n, m £ N. This process can 
be described equivalently as a time-changed M/M/l queue (see Proposition [3J]) or 
as a sequence of independent branching processes (see Proposition [321). The time- 
change is the discrete analog of the Lamperti transform between continuous-state 
branching processes and Levy processes, see Lamperti [8]. As it will be seen these 
two viewpoints are complementary. 

Let {<J n ) be the sequence of birth instants (i.e., positive jumps) of (Z(t)) and 
(B a (t)) the corresponding counting process of (<J n ), for t > 0, 



i>i 



Proposition 3.1 (Queueing Representation). If Z(0) = zetf, then 

(5) (Z(t),t>0) d ^ L (L(C(t)),t>0), 

where (L(t)) is the process of the number of jobs of an M/M/l queue with input 
rate \xz and service rate v and with L(0) = z and C(t) = inf {s > : A(s) > t}, 
where 

r* i 

A(t) = / — - du. 

y> J 1VL(«) 

Proof. It is not difficult to check that the process (M(t)) d =' (L(C(t))) has the 
Markov property. Let Qm be its Q-matrix. For z > 0, 

¥(L(C(h)) =z + l\ L(0) = z)= n z E(C(h)) + o{h) = n z {z V l)h + o(h), 

hence qniz, z + 1) = ^z( z V 1). Similarly (?m( z j z — 1) = vz. The proposition is 
proved. □ 

Corollary 3.1. For any 7 > (fiz - v) V and z = Z(0) e N, 

(6) E z (j2e~-« a *^J < +00. 

Proof. Proposition 13.11 shows that, in particular, the sequences of positive jumps 
of (Z(t)) and of (L(C(t))) have the same distribution. Hence, if Ma Z — (tn) is the 
arrival process of the M/M/l queue, a Poisson process with parameter fiz, then, 
with the notations of the above proposition, the relation 

K) (A(t n )) 
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holds. By using standard martingale properties of stochastic integrals with respect 
to Poisson processes, see Rogers and Williams [IS] , one gets for t > 0, 

= E Z QT e-^M^ds)^ =fi z E z Ijf^e^d. 

/>OC 

(7) = /iz / e- T "E, (Z(u) V 1) du, 

Jo 

where Relation §5§ has been used for the last equality. Kolmogorov's equation for 
the process (Z(t)) gives that 

<P(t) d = E z (Z(t)) =fjtz / E z V 1) - v / E z (Z(tt)) du 



< (nz — v) / du + /izi, 

JO 

therefore, by Gronwall's Lemma, 

(j){t) < 0(0) + nz / ue^*-"> tt dit < z + /iz te^*~^*. 
Jo Hz -v 

From Equation ([7]), one concludes that 

e z (x>~ 7<T ") = E - (x>~^ (tn) l < +°°- 

The proposition is proved. □ 

Before hitting 0, (Z(t)) can be seen a Bellman-Harris branching process with 
Malthusian parameter a = Hz — v, see Athreya and Ney [3]. This Bellman- Harris 
branching process describes the evolution of a population of independent particles: 

def 

each particle, at rate A =' Hz + f, cither splits into two particles with probability 
P *== Hz I {.Hz + v) or dies with probability 1 — p. These processes will be referred 
to as (p, A)-branching processes in the sequel. 

A (p, A)-branching process survives with positive probability only when p > 1/2, 
in which case the probability of extinction q is equal to q = (1 — f>)/f> = vj Hz- The 
process (Z(t)) only differs from a (p, A)-branching process insofar it regenerates 
after hitting 0, after a time exponentially distributed. When it regenerates, it 
again behaves as a (p, A)-branching process (started with one particle), until it 
hits again. 

Proposition 3.2 (Branching Representation). If Z(0) = z G N and {Z(t)) is a 
(p, X) -branching process started with z£N particles and T its extinction time, then 

(Z(t),0 <t<T) dl = (Z(t),0 <t<f), 

where T = mi{t > : Z(t) = 0} is the hitting time o/O by (Z(t)). 

Corollary 3.2. Suppose that hz > v. Then ¥ z -almost surely for any z > 0, there 
exists a finite random variable Z[po) such that, 

lim e-^ z - v)t Z(t) = Z(oo) and Z(oo) > 0. 
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Proof. When fiz > v, the process {Z(t)) couples in finite time with a supercritical 
(p, A)-branching process (Z(t)) conditioned on non-extinction; this follows readily 
from Proposition 13.21 (or see the Appendix for details). Since for any supercritical 
(p, A)-branching process, (exp(— (fiz — v)t)Z{t)) converges almost surely to a finite 
random variable Z(oo), positive on the event of non-extinction (see Nerman 
one gets the desired result. □ 

Due to its technicality, the proof of the following result is postponed to the 
Appendix; this result is used in the proof of Proposition [ 

Proposition 3.3. Suppose that jiz > v. If 



2 - x - JxU - 3x) 
(8) r,*(x) = ^0<*<1, 

then for any < rj < r\*(yj 'nz), 



sup 

z>0 



e z (sup (e^-^B^tyvyj 



< 



A Yule Process Killed at the Birth Instants of a Renewing Bellman- 
Harris Process. We now consider the main interacting branching processes: in 
addition to (Z(t)), one considers an independent Yule process (Y(t)) with parameter 
fiw (its Q-matrix is defined by Relation ((T|) with fi = fiw). We proceed similarly as 
in Section[2] a process (W(t)) is defined by killing one individual of (Y(t)) at each 
of the birth instants (<r„) of (Z(t)) (see Alsmeyer [2] and the references therein for 
related models). Recall that (W(t)) is given by Formula <j2j> . The following results 
are key to analyzing the two-chunk network of Section [5l 

Proposition 3.4. Assume that fiz — v > fJ-w > an d let Hq be the extinction time of 

(inn), i.e., 

Ho = inf{t > : W(t) = 0}, 
then the random variable Hq is almost surely finite and: 

(i) Z(H ) - Z(0) < e^ vHo M^ where 

Mp = sup (e-^ wt Y(t)) . 
t>o 

(ii) There exists a finite constant C such that for any z > and w > 1, 

(9) E (W!Z) (H )<C(log(w) + l). 

In ^ the subscript (w, z) refers to the initial state of the Markov process 
(W(t), Z(t)). More generally in the rest of the paper we use the convention that if 
(U(t)) is a Markov process then the index u of P u and E u will refer to the initial 
condition of this Markov process. 

Proof. Define a = pz — v. Concerning the almost sure finiteness of Hq, note that 
Equation © entails that W(t) < Y(t)-B a (t) for all t > on the event {H = +oo}. 
As t goes to infinity, both exp(— /j,wt)Y(t) and exp(— at)B a (t) converge almost 
surely to positive and finite random variables (see Nerman which implies, 

when a — [iz — v > pw, that W(t) converges to — oo on {H = +oo}, and so this 
event is necessarily of probability zero. 

The first point (i) of the proposition comes from Identity © evaluated at t = H : 

(10) Z(H ) - Z(0) < B a (H ) < Y(H ) < e^ vH °Mp. 
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By using the relation exp(a;) > x, Equation ([9]) follows from the following bound: 
for any 77 < rf{yj[iz) (recall that rf is given by Equation ((HJ), 



(11) sup L-"E Kz) ( e *f*-*w)B \ 

w>l,z>0 L \ / 



< +OO. 



So all is left to prove is this bound. Under P( ti , j2 ), (Y(t)) can be represented as the 
sum of w i.i.d. Yule processes, and so M Y < My 1 + • • • + M y ^ with (M Y A i.i.d. 
distributed like My under P(i, z ); Inequality (fTCTI) then entails that 

e (a- w )ffo < [ W yj ) x sup (e at /B ff (t)) . 
By independence of (MyJ and (B a (t)), Jensen's inequality gives for any 77 < 1 

hence the bound (fTTj) follows from Proposition 13.31 □ 

One concludes this section with a Markov chain which will be used in Section [5] 
Define recursively the sequence (V n ) by, Vq = v and 

A n {v n ) 

(12) V n+1 = T n,k, n>Q, 

k=l 

where for each n (I n _k ,k>l) are identically distributed integer valued random vari- 
ables independent of V n and A n (V n ), and such that E(/„,i) = p for some p £ (0, 1). 
For v > 0, A n (v) is an independent random variable with the same distribution as 
Z(H ) under P(i,„), i.e., with the initial condition (W(0),Z(0)) = (l,u). 

The above equation (|12p can be interpreted as a branching process with immi- 
gration, see Seneta [17] . or also as an auto- regressive model. 

Proposition 3.5. Under the condition fiz — v > [iw , if (Yn) * s ihe Markov chain 
defined by Equation (|12p and, for K > 0, 

N K = mi{n > : V n < K}, 

then there exist 7 > and K £ N such that 

(13) E(N K \V = v) < -log(l + v), Vtj>0. 

7 

The Markov chain (V n ) is in particular positive recurrent. 

Proof. For Vq = v £ N, Jensen's Inequality and Definition (fT2")l give the relation 



(14) ( ^-y) 



1 + pZ(H ) 
1 + v 

From Proposition 13 .41 and by using the same notations, one gets that, under P(i jlI ), 

Z{H Q ) < v + e» wHo M£, 

where (Y(t)) is a Yule process starting with one individual. By looking at the birth 
instants of (Z(t)), it is easily checked that the random variable Hq under P(i jt ,) is 
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stochastically bounded by Hq under P(i,o)- The integrability of Hq under P(i,o) 
(proved in Proposition I3.4[) and of My give that the expression 

, 'l+p(v + ef" nHo Mt) 
lor 1 



1 + V 

bounding the right hand side of Relation (|14l) is also an integrable random variable 
under P(i,o)- Lebesgue's Theorem gives therefore that 

1 + Vi ' 



lim sup 

v— >+oo 



log 



1 



< \ogp < 0. 



Consequently, one concludes that v i— > log(l + v) is a Lyapunov function for the 
Markov chain (V n ), i.e., if 7 = — (logp)/2, there exists K such that for v > K, 

E v log (1 + Vi) - log (1 + v) < -7. 

Foster's criterion, see Theorem 8.6 of Robert [15], implies that (V n ) is indeed ergodic 
and that Relation (fT3| holds. □ 

4. The Single-Chunk Network 

This section is devoted to the study of a class of two-dimensional Markov jump 
processes (Xo(t), Xi(t)), the corresponding Q- matrix Q r is given, for x = (x<j, x\) £ 
N 2 , by 

!Vt r [{x Q ,x\), (x + l,xi)] = A, 

tl r [(xo,xi), {x - l,xi + 1)] = /ir(io,ai)(a:i V l)l{a; D >o}, 
n r [(x ,a;i), (x ,a;i - 1)] = ^i, 

where x 1— > r(x), referred to as the rate function, is some fixed function on N 2 with 
values in [0,1]. 

From a modeling perspective, this Markov process describes the following system. 
Requests for a single file arrive with rate A, the first component Xo(t) is the number 
of requests which did not get the file, whereas the second component is the number 
of requests having the file and acting as servers until they leave the file-sharing 
network. The constant /1 can be viewed as the file transmission rate, and v as 
the rate at which servers having the file leave. The term r(xo,xi) describes the 
interaction of downloaders and uploaders in the system. The term x\ V 1 can be 
interpreted so that there is one permanent server in the network, which is contacted 
if there are no other uploader nodes in the system. A related system where there 
is always one permanent server for the file can be modeled by replacing the term 
X\ V 1 by X\ + 1. See the remark at the end of this section. 

Several related examples of this class of models have been recently investigated. 
The case 

r(x ,xi) = ■ 

X + X\ 

is considered in Nuhez-Queija and Prabhu [13] and Massoulie and Vojnovic [10] : 
in this case the downloading time of the file is neglected. Susitaival et al. [19] 
analyzes the rate function r{x) 

1 \ 1 f x o 

r(xo,xi) = 1 A a — 

V x i 
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with a > and a A b denotes min(a, b) for a, b € M. This model allows to take into 
account that a request cannot be served by more than one server. See also Qiu and 
Srikant Q3]. 

With a slight abuse of notation, for < S < 1, the matrix fl$ will refer to 
the case when the function r is identically equal to S. Note that the boundary 
condition x\ V 1 for departures from the first queue prevents the second coordinate 
from ending up in the absorbing state 0. Other possibilities are discussed at the 
end of this section. In the following {X r {t)) = (X£{t),X[(t)) [resp. (X s (t))] will 
denote a Markov process with Q-matrix f2 r [resp. £1$]. 

Free Process. For 5 > 0, Qs denotes the following Q-matrix 

(Qs[(y,z),(y + l,z)} = A, 

(16) (Q 5 [(y,z),(y-l,z + l)] =*V1), 
[Qs[{y,z),(y,z- 1)] = vz. 

The process (Y s (t)) = (Y^(t),Yf(t)), referred to as the free process, will denote 
a Markov process with Q-matrix Q$. Note that the first coordinate (Yg(t)) may 
become negative. The second coordinate (Yf(t)) of the free process is the birth- 
and-death process with the same distribution as the process (Z(t)) introduced in 
Section |3l It is easily checked that if ps defined as 5p/v is such that ps < 1, then 
(Yi(t)) is an ergodic Markov process converging in distribution to Yi(oo) and that 

(17) \*(6) i/E(y*(oo)) = ^E(y/(oo) V 1) = ^ -. 

(l-p s ){l - log(l - ps)) 

When ps > 1, then the process (Yi(i)) converges almost surely and in expectation 
to infinity. In the sequel A*(l) is simply denoted A*. 

In the following it will be assumed, see Condition (C) below, that the rate function r 
converges to 1 as the first coordinate goes to infinity; as will be seen, the special 
case r = 1 then plays a special role, and so before analyzing the stability properties 
of (X r (t)), one begins with an informal discussion when the rate function r is 
identically equal to 1. Since the departure rate from the system is proportional to 
the number of requests/servers in the second queue, a large number of servers in 
the second queue gives a high departure rate, irrespectively of the state of the first 
queue. The input rate of new requests being constant, the real bottleneck with 
respect to stability is therefore when the first queue is large. The interaction of the 
two processes (X^t)) and (Xl(t)) is expressed through the indicator function of 
the set {Xo(t) > 0}. The second queue (X[(t)) locally behaves like the birth-and- 
death process (Yi(t)) as long as (Xg(i)) is away from 0. The two cases pi > 1 and 
pi < 1 are considered. 

If p\ > 1, i.e., p > v, the process (Xj(i)) is a transient process as long as 
the first coordinate is non-zero. Consequently, departures from the second queue 
occur faster and faster. Since, on the other hand, arrivals occur at a steady rate, 
departures eventually outpace arrivals. The fact that the second queue grows when 
(Xo(t)) is away from stabilizes the system independently of the value of A, and 
so the system should be stable for any A > 0. 

If pi < 1, and as long as (X (t)) is away from 0, the coordinate (Xl(t)) locally 
behaves like the ergodic Markov process (Y±(t)). Hence if (Xg(i)) is non-zero for 
long enough, the requests in the first queue see in average E(Y 1 1 ((>o) V 1) servers 



INTERACTING BRANCHING PROCESSES 



11 



which work at rate fi. Therefore, the stability condition for the first queue should 
be 

A < (po) V 1) = A* 

where A* = A*(l) is defined by Equation (fl7|) . Otherwise if A > A*, the system 
should be unstable. 

Transience and Recurrence Criteria for (X r (t)). 

Proposition 4.1 (Coupling). If X r (0) — F 1 (0) G N 2 , there exists a coupling of 
the processes (X r (t)) and (Y (t)) such that the relation 

(18) X r Q (t) > Y^t) and X[(t) < *?(«), 

holds for all t > and for any sample path. 
For any < 5 < 1 , if 

t s = inf{t > : r(X r {t)) < 6} and a = inf{t > : X£(t) = 0}, 

and if X r (0) — Y s (0) £ N 2 then there exists a coupling of the processes (X r (t)) 
ana (Y s (t)) such that, for any sample path, the relation 

(19) Xj(t) < Y S (t) and X{{t) > F/(t) 
holds for all t < t$ A a . 

Proof. Let X r (0) — (xq,Xi) and y 1 (0) = (y, z) be such that xq > y and X\ < z, 
one has to prove that the processes (X r (t)) and can be constructed such 

that Relation (IT8l) holds at the time of the next jump of one of them. See Leskela 
for the existence of couplings using analytical techniques. 

The arrival rates in the first queue are the same for both processes. If x\ < 
z, a departure from the second queue for (Y" 1 (t)) or (X r (t)) preserves the order 
relation (|18[) and if x\ = z, this departure occurs at the same rate for both processes 
and thus the corresponding instant can be chosen at the same (exponential) time. 
For the departures from the first to the second queue, the departure rate for (X r (t)) 
is fir(xo, x\){x\ Vl)l{ :Eo> o} < (i(zVl) which is the departure rate for (F 1 (t)), hence 
the corresponding departure instants can be taken in the reverse order so that 
Relation (|I8j) also holds at the next jump instant. The first part of the proposition 
is proved. 

The rest of the proof is done in a similar way: The initial states X r (0) — (xq, x\) 
and Y s (0) = (y,z) are such that xq < y and x\ > z. With the killing of the 
processes at time t$ A a one can assume additionally that xo ^ and that the 
relation r(xo,Xi) > 5 holds; Under these assumptions one can check by inspecting 
the next transition that (fT9"| holds. The proposition is proved. □ 

Proposition 4.2. Under the condition /j, < v, the relation 

liminfSl > A -A* 

t^+oo t 

holds almost surely. In particular, if fi < v and X > X* , then the process (X r (t)) is 
transient. 

Proof. By Proposition 14. R one can assume that there exists a version of (F 1 (t)) 
such that X£(0) = 5^(0) and the relation > Y^(t) holds for any t > 0. From 

Definition (|T6| of the Q-matrix of (F 1 (t)), one has, for t > 0, 

Y o 1 (t) = Y a 1 (0)+M x (t)-A(t), 
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where (Af\(t)) is a Poisson process with parameter A and (A(i)) is the number of 
arrivals (jumps of size 1) for the second coordinate (Yj (t)): in particular 



E(A(t)) = fiE (J^ Y^(s) V 1 ds^ 



Since (Yi(i)) is an ergodic Markov process under the condition \i < is, the ergodic 
theorem in this setting gives that 

lim -A(t) = lim -E(A(t)) = fiE (Y}(oo) V l) = A*, 

t—>-\-oo t t—^-\-oo t 

by Equation (IF]), hence (Yj(t)/t) converges almost surely to A — A*. The proposi- 
tion is proved. □ 

The next result establishes the ergodicity result of this section. 
Proposition 4.3. If the rate function r is such that, for any x± G N, 
(C) lim r(x ,xi) = 1, 

Xq— > + OC 

and if fj, > v, or if fi < is and A < A* with 



(20) A *-(l-p)(l-log(l-p))' 
and p = (a/v, then (X r (t)) is an ergodic Markov process. 

Note that Condition (C) is satisfied for the functions r considered in the models 
considered by Nufiez-Queija and Prabhu [13] and in Susitaival et al. [19.. See 
above. 

Proof. If x — (xq,xi) g M 2 , \x\ denotes the norm of x, \x\ = \xq\ + \xx\. The 
proof uses Foster's criterion as stated in Robert [T5] Theorem 9.7]. If there exist 
constants Kq, Ki, to and t\ such that 

(21) sup E (x0)Xl) (|X r (t 1 )|-|a;|) < 0, 

n>Xi 

(22) sup M (xOi:Cl) (|^(£ )| - < 0, 

x 1 <K 1 ,x 2 >K 2 

then the Markov process (X r (t)) is ergodic. 

Relation ([21]) is straightforward to establish: if xi > Kx, one gets, by considering 
only K\ of the x\ initial servers in the second queue and the Poisson arrivals, that 

E (x0jXl) (\X r (l)\-\x\)<X-K 1 (l-e- v ), 

hence it is enough to take t\ = 1 and K\ = (A + 1)/(1 — e _ly ) to have Relation ([2~T]) . 

One has therefore to establish Inequality ([22]) . Let t$ and a be the stopping 
times introduced in Proposition |4T] one first proves an intermediate result: for any 
t > and any x\ G N, 

(23) lim F (xOtXl) (aAT 5 <t)=0. 

Xq— > + 0O v ' 

Fix x\ G N and t > 0: for e > 0, there exists D\ such that 

sup Yl{s) >Di) < £, 

0<s<* / 
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from Proposition l4.il this gives the relation valid for all xq > 0, 

P(*o,*0 ( SU P ^ D ^ £ - 

\0<s<t / 

By Condition (C), there exists 7 > (that depends on x\) such that r(xo,xi) > 5 
when xo > 7. As long as (X r (t)) stays in the subset {(yo,Vi) ■ Ui < -Di}, the tran- 
sition rates of the first component (XJ(i)) are uniformly bounded. Consequently, 
there exists K such that, for xq > K, 



supXJ(s) < 7, supX[(s) < D u 

s<t s<t 



< 



Relation ([25)1 follows from the last two inequalities and the identity 

P(*o,*i)( ff A r 5 < <) < ¥ (X0}X1) (fpXo( s ) < 7j • 

One returns to the proof of Inequality ([212]) . By definition of the Q-matrix of the 
process (X r (t)), 



E (x0iXl) (\X r (t\) - |a:|) =Xt-v E (x0iXl) (X{(u))du, ieN , t >0. 

Jo 

For any x £ N 2 , there exists a version of (Y 5 {t)) with initial condition Y 5 (0) 
X r (0) = x, and such that Relation (|T5)) holds for i < ra A a, in particular 

E K (X[(t)) >E x (Xr(t);t<r 5 Aa) 

> E x (y/(t);f <tj Ait)= (!?(*)) - E x (F/(t);t > r 5 A a). 

Cauchy-Schwarz inequality shows that for any t > and i£N 2 



^(^/(w); r 5 A a < u) du < 



< y/Vjji Act < t) 



(y'(tt))' 



by gathering these inequalities, and by using the fact that the process {Yf(t)) 
depends only on x\ and not xq, one finally gets the relation 

1 



(24) -W. x {\X{i)\-\x\)<\-- I E X1 (Yf(u)) du + c( Xl ,t)y/¥ x (Ts A a < t) 



with 



c(xi,t) 



t 



E, 



(Yf(u)Y 



du. 



Two cases are considered. 

(1) If fj, > u, if S < 1 is such that > v, the process (F* (t)) is transient, so 
that 



lim — 

t— S- + OC £ 



E^F/^cfe = +00, 



for each x\ > 0. 
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(2) If /i < v, one takes 6 = 1, or if \i = v, one takes 5 < 1 close enough to 1 so 
that A < A* (6). In both cases, A < A*(<5) and the process (Y-f(t)) converges 
in distribution, hence 

pi 



lim - / E xl (Yf(u))du = vE (Yf(oo)) = X*(S) > A 
t^+oc t J 

for each x\ > 0. 

Consequently in both cases, there exist constants r\ > 0, S < 1 and to > such that 
for any X\ < K\, 

(25) A-f (\ Xl {Y*{u))du<-r,, 

to Jo 

with Relation ([24| . one gets that if x\ < K\ then 



^E x (\X(t )\ - |ar|) < -r? + c*V>*(t4 A <r < i ), 
to 

where c* = max(c(n, to), < n < K%). By Identity ([23|l . there exists i^o such that, 
for all xo > Kq and xi < ifi, the relation 



y^Wofa A<7<* ) 



<5 

- 2 



holds. This relation and the inequalities (|2"5"]) and ([M]) give Inequality (|22[) . The 
proposition is proved. □ 

Another Boundary Condition. The boundary condition xi V 1 in the transition 
rates of (JC(t)), Equation (|15[) . prevents the second coordinate from ending up in 
the absorbing state 0. It amounts to suppose that a permanent server gets activated 
when no node may offer the file. Another way to avoid this absorbing state is to 
suppose that a permanent node is always active, which gives transition rates with 
xi + 1 instead. This choice was for instance made in Nunez-Queija and Prabhu [13) . 
All our results apply for this other boundary condition: the only difference that is 
when v > /x, the value of the threshold A* of Equation is given by the quantity 
A* = [iv I (y — fj). 



5. The Two-Chunk Network 

In this section it is assumed that a file of two chunks is distributed by the file- 
sharing system corresponding to Figure [TJ Chunks are delivered in the sequential 
order. The analysis of this simple file-sharing system gives a hint of the difficulties 
one can encounter when dealing with multiple chunks, in particular with nodes 
distributing and receiving chunks. 

For k = 0,1 and t > 0, the variable Xj~{t) denotes the number of requests 
downloading the (fc+l)st chunk; for k = 2, A~2(i) is the number of requests having 
all the chunks. When taking into account the boundaries in the transition rates 
described in Figure HJ one gets the following Q-matrix for the three-dimensional 
Markov process (X(t)) = (X k (t), k = 0, 1, 2) 

Q(f)(x) = X[f(x + e ) - /Or)] + m{x x V l)[f(x + e t - e ) - f(x)]l {xo>Q} 

+ fi 2 (x2 V l)[f(x + e 2 - ei) - f(x)]l {xi>0} + vx 2 [f{x - e 2 ) - f(x)], 
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where x 6 N 3 , / : N 3 ->■ E+ is a function and for, < k < 2, eu € N 3 is the kth unit 
vector. Note that, as before, to avoid absorbing states, it is assumed that there is 
a server for the kth chunk when Xk = 0. 

The stability behavior of (X(t)) depends on the values of the parameters fii, /i 2 
and v. Three cases need to be distinguished: (1) /ii > /i 2 — v > 0, (2) fi 2 — v > fi\ 
and (3) /z 2 < v. In each case, the method of proof is similar to the one used in 
the proof of Proposition 14.31 and relies on Foster's criterion as stated in Robert [H2 
Theorem 9.7]. One investigates the case when the size xq + x\ + x 2 of the initial 
state X(0) = x = {xq,X\,X2) is large. We briefly discuss the three situations which 
will be analyzed in the following. 

(1) If X2 is large, then the total number of customers decreases instantaneously. 

(2) If x\ is large and X2 is small, then the last queue (X 2 (t)) behaves like the 
birth-and-death process (Z(i)) of Section [3l which is transient if [12 > v 
and stable otherwise. 

(3) If x is large and x\ and x 2 are small, the two last queues (Xi(t), X 2 (t)) 
behave like the system (X s (tj) — [X^{t), X^t)), defined by its generator 



(26) Q s (f)(z) = ^(z! V l)[f(z + ei) - /(*)] 

+ /i 2 (z 2 V l)[/(z + e 2 - ei) - f(z)]l {zi>0} + vz 2 [f{z - e 2 ) - /(*)], 



where z £ N 2 and / : N 2 — > K+ is an arbitrary function. 

Since (X s (£)) corresponds to the two last queues when the first queue is satu- 
rated, we call it the saturated system: it plays a crucial for the stability of (X(t)), 
in particular the asymptotic behavior of E(X 2 s (t)) determines the output of the 
system when the first queue is saturated. The properties of (X s (t)) that will be 
needed are gathered in the next proposition. 

Proposition 5.1 (Saturated system). Let x €N 2 . 

(1) If fit > [i 2 - v > 0, then E x (X$(t)) -> +00 as t -}■ +00. 

(2) If I-I2 — v > fii, then (X s (t)) is positive recurrent. 

(3) If v > fj, 2) then ^E x (X^(t)) — > A* as t — > +00, with X* given by 



Proof. The crucial observation is that as long as Xf(t) > 0, the process (X s (t)) 
can be coupled with the process (W(t), Z(t)) of the last part of Section [3l i.e., 
(Xf (i)) behaves like a birth-and-death process and (Xf (t)) is a Yule process with 
particles killed at the instants of birth of (X^(t)). 

More formally, in the sequel let (Z(t)) be the process with Q-matrix defined by 
Relation (0} with fiz = M2 5 an d define (er„) its sequence of birth times. Let (W(t)) 
be a Yule process with parameter [i\ killed at (<r n ). Then (X s (t)) and (W(t), Z{t)) 
can be coupled so that 



with T — inf{i > : Xf{t) — 0}. We give a full proof of the two first cases, the 
last case uses similar techniques and so we leave the details to the reader. 

Proof of case (1). Since /i 2 > u, the process (e^^^^Z^)) converges almost surely 
to a finite and positive random variable Mz (00) by Corollary [321 Moreover, since 



(27) 



A* = 



(l-Ate/fXl-logtl-Aia/f))' 



(t), < t < T) = (W(t), Z(t), < t < T) 
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Mi > A*2 — v > 0, Corollary 13. II entails that 

g-Mi^n < ^qq^ almost surely. 

n>l 

By Proposition 12.11 this shows that (W(t)) survives with positive probability. In 
this event, say £, it holds that W(t) > 1 at all times and 

lim Z(t) — +oo. 

In the event £, W(t) > 1 and therefore (Xf(i), Xf (i)) = (W(t), Z(t)): in particu- 
lar, we have exhibited an event of positive probability where X^t) — > +oo, which 
proves the claim. Note that we implicitly assumed that no coordinate of the initial 
state is zero, these cases being dealt with easily. 

Proof of case (2). This is the most delicate case. By Proposition 13.41 and since 
[12 — v > /ii, one has that (Xf(t)) returns infinitely often to 0. When (Xf(t)) is 
at it jumps to 1 after an exponential time with parameter hi. One defines the 
sequences (i?o,n)> (Efi,n) an d (S n ) by induction: So = by convention, and for 
n > I, S n = J2k=i(Ho,k + E^.k), 

H ,n = inf{< > S n -! : Xf (t) - 0} and E^ n = M{t > S„_i + H 0}71 : Xf (t) = 1}. 

The sequence (E^n) represents the successive times at 0, (S n ) the times at which 
(Xf (t)) hits 0, and (-ffo.n) the times needed for (Xf (i)) to hit 0. On the time inter- 
vals [S n ,S„ + iJo,n+i]) {X s (t)) can be coupled with a version of {W(t), Z(t)). This 
shows that, with the notations of Proposition 13.41 Xf (SVi) is equal in distribution 
to Z(H ) with the initial condition (W(0),Z(0)) = (Xf (5„_i), Xf^n-i)). Then 
while Xf is at 0, the dynamics of Xf is simple, since it just empties. These two 
remarks show that X^Sn+i) is equal in distribution to 

/Z(H ) 

(Xi(S n+1 )\X s (S n )=x) d ^ I J2 t { E^>E n }\(W(0),Z(0)) = x 

with (E u j) i.i.d., exponentially distributed with parameter v and independent of 
E^. For i > 1, ii 1 ,,^ represents the residence time of the ith customer present in 
the second queue at time S n + Ho, n +i- 

Consequently, and since Xf(5„) = 1 for n > 1, (Xf (S n ),n > 1) is a Markov 
chain whose transitions arc defined by Relation (1121) with p = u / (u + fi\); note that 
(X^iSn), n > 0) has the same dynamics only when Xf (0) = 1. In the sequel define 
for any K > the stopping Tk- 

T K = inf{n > : X§(S n ) < K,Xf(S n ) = 1}. 

The ergodicity of (X s (t)) will follow from the finiteness of E/ X1}X2 -\(Tk) for some K 
large enough and for arbitrary x — (x-y,X2) € N 2 : indeed, Tk is greater than the 
time needed to return to the finite set {x\,X2 ■ X\ = 1,X2 < X}. The strong 
Markov property of (X s (i)) applied at time Si gives 

and so one only needs to study T K conditioned on {Xf (0) = 1} since E( Xl ^ 2 )(S'i) 
is finite in view of Proposition 13.41 Thus one can assume that (Xf (S n ),n > 0) is 



^( Xi ,x2)( t k) < 2E( Xlia . 2 )(Si) +E (a , liX2 ) 
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a time- homogeneous Markov chain whose transitions arc defined by Relation (|12[) . 
and with Nk defined in Proposition 13.51 the identity 

N K 

(28) T K =Y J {Ho, i + E^ i ) 

holds. For i > 0, the Markov property of (X% (S n ),n > 0) gives for any x 2 > 

E (l,x 2 ) (Ho,it{i<N K }) = E (l,z 2 ) ( E (l,X 2 s (Si)) (#o) l{i<iV K }J ■ 

With the same argument as in the proof of Proposition 13.51 one has 

E (i,xi( St ))( H o) < %,o)(#o) < +oo, 

with Equations (|28p and (|TT?)) of Proposition 13.51 one gets that for some 7 > and 
some K > 0, 

E (XUX2) (T K ) < 2E {XUX2) (S 1 ) + CE {XUX2) [log(l + X 2 s (^))] 

with the constant C = (En ) o)(i?o)+l/j u 2)/7- This last term is finite for any (x±, 
in view of Proposition 13.41 which proves the proposition. 

Sketch of proof of case (3). Since fii > > f/,2 ~ v . Corollary 13.11 shows that 
E n >i e "'"°" 1S nn he almost surely, and Proposition 12.11 implies that (Xf(£)) sur- 
vives with positive probability. Using the fact that the process (Z(t)) is stable since 
/i2 < v, together with similar arguments as before, one can show that the process 
(Xf(t)) actually survives almost surely. Thus after some time {X^it)) can be cou- 
pled forever with (Z(t)), and A* is precisely the value of i/E(Z(oo)) in stationary 
regime. □ 

Proposition 5.2. Let (X(t)) = (Xo(t) ) Xi(t),X2(t)) be the Markov process de- 
scribing the linear file- sharing network with parameters A,/zi,/^2 and v. 

(1) If fii > fj,2 — v > 0, then (X(t)) is ergodic for any A > 0. 

(2) If [12 —v>fJ,i, let X s be the value of the expectation of vX 2 at equilibrium, 
which exists by Proposition 15. 1 1 Then (X(t)) is ergodic for A < X s and 
transient for X > X s . 

(3) If v > /12, then (X(t)) is ergodic for X < X* and transient for X > X* , with 
X* given by ([27]). 

In the second case, preliminary investigations seem to suggest that X s = +00 
is possible. In such situations, the stability region would actually be infinite. The 
intuitive reason is that (Xg (£)) grows exponentially fast in periods where Xf (t) > 0, 
and for a duration H . Thus when (Xf(t)) hits 0, X§(H Q ) w e^ 2 '^" , and 
although some exponential moments of H are finite by (|11[) . it seems plausible 
that E(e^ 2 ~^ ffo ) = +00 sometimes. 

Proof. In all three cases the proof relies on Foster's criterion as stated in Theo- 
rem 9.7 of Robert [15]. Let X(0) = x = (x ,xi,x 2 ) G N 3 , assume that = 
xq + x\ + X2 is large. First we inspect the case when x 2 is large, it is common to 
all the cases. 

Since the last queue serves at rate v each request, for t > 0, 

Ex(||a:(*)||) < IMI + At - X2 (1 - e-^) . 
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Define t-x = 1 and let K2 be such that A — if 2(1 — e ") < — 1, then the relation 

E x (]|x(t 2 )||)- W<-i, 

holds when X2 > K2, for all xo,xi > 0. In the sequel the constants £2 and if 2 are 
fixed, and the process (Z(t)) is defined similarly as in the proof of Proposition 15. II 
We give a full proof of case (1), the other cases use similar techniques and so we 
leave the details to the reader. 

Proof of case (1). Assume that /ii > fi2 — v > 0: in particular (Z(t)) is transient, 
thus there exists some ti such that for any X2 < if 2, 

v I E X2 (Z(u)) du > A^ +2. 
Jo 

The two processes (Z(t)) and (X(t)) can be built on the same probability space 
such that if they start from the same initial state, then the two processes (Z(t)) 
and (AT2(t)) are identical as long as X\{t) stays positive. Since moreover the hitting 
time inf{t > : X\(t) = 0} goes to infinity as x\ goes to infinity for any X2 < K%, 
one gets that there exists K\ such that if x\ > K\ and X2 < K2, then the relation 

EsfllXXti)!!) - |M| = Xh - v [ 1 E x {X 2 {u))du 

Jo 

<\ti-(vj^ E X2 {Z(u)) du - l\ < -1 

holds. 

In a similar way, by saturating the first queue and coupling the two last queues 
(Xi(t), X2{t)) with the saturated system (X s (tj), one gets in view of case (1) of 
Proposition 15.11 that there exist constants to and Kq such that if X2 < K2 , and 
x\ < K\ and xq > Kq, then 

E x (||X(t )||)-||x|| <-l, 
which ends the proof of this case. 

Sketch of proof of case (2). Since \i2 > v, the process (Z(t)) is still transient and 
one can choose the same constants t\,K\ as in the previous case to deal with the 
case x\ > K\,X2 < K2- When xq > Kq and x\ < K\, X2 < ^2, then one can justify 
for Kq and to large enough the approximation 

E x {\\X{t )\\) - Hzll « At - v f ° E (xuX2) (Xi(u)) du^(X- X s )t . 

Jo 

Hence for A < X s , the stability result follows. For A > A , this suggests that the 
first coordinated (Xo(t)) builds up, and similarly as in the proof of Proposition 14.21 
one can show that indeed in this case 

X (t) 

liminf > A — X s > 0, almost surely. 

*->+oo t 

Sketch of proof of case (3). Since v > fi2, the process (Z(t)) is stable, and it acts as 
a bottleneck on the system. When Xo < Kq and x\ > K±, then one can as before 
justify the approximation, for K\ and t\ large enough, 

E^WX^W) - \\x\\ « Atr - v f 1 E X2 (Z(u)) du^(X- A*)ti. 

Jo 
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Similarly, the third case of Proposition 15.11 shows that when xq > Kq and x\ < 
Ki,X2 < K2, again for Kq and to large enough, 

f*0 



E x (\\X(t )\\)-\\x\\n\t - V / E (xuX2) (XZ(u))du*(\-\*)t . 

Jo 

This shows that the system is ergodic for A < A*. The transience in the other case 
A > A* follows easily from the fact that the last queue (^(t)) is stochastically 
dominated by (Z(i)), thus A* is the maximal output rate of the system. □ 



Appendix A. Proof of Proposition 13.3 



In this appendix the notations of Section [3J are used. Since the random variable 
(B a (t) I Z(0) = 0) is stochastically smaller than (B a (t) \ Z(0) = z) for any z E N, 
it is enough to show that for r\ < rfiyj ^z) 



En 



sup {f at B a {t)^) 

,t>CTl 



< +OO, 



where a = [iz — v > 0. 

Note that the process (B a (t + ai),t > 0) under F has the same distribution as 
{B a {i) + l,t > 0) under Pi, and by independence of <j\, an exponentially random 
variable with parameter /iz, and (B a (t + a\),t > 0), one gets 



sup {e r > at B a {t)- r >) 

.t>ai 



= E (e r < 



IE, 



sup 

t>0 



(e" Qt (B^t) 



1) 



Since a < fiz and r]*(v//j,z) < 1, then Eo (exp^acri)) is finite, and all one needs 
to prove is that the second term is finite as well. 
Define r as the last time Z(t) = 0: 

r = sup{i > : Z(t) = 0}, 

with the convention that r = +00 if (Z(t)) never returns to 0. Recall that, because 
of the assumption [iz > with probability 1, the process (Z(t)) returns to a 
finite number of times. 

Conditioned on the event {t = +00}, the process [Z(t)) is a (p, A)-branching 
process conditioned on survival, with A = fiz + v and p = /j,z/X. Such a branching 
process conditioned on survival can be decomposed as Z — Zru + Y, where (Y(t)) 
is a Yule process (Y(t)) with parameter a. See Athreya and Ney [3]. Consequently, 
for any < i] < 1, 



Ei 



sup 

t> 



P ( e " Qt (B a (t) + !)-*>) 
v ' 



< Ei 



sup (e" Qt r(t)-") 
t>o 



Since the nth split time t n of (Y(t)) is distributed like the maximum of n i.i.d. 
exponential random variables, Y(t) for t > is geometrically distributed with 
parameter 1 — e~ Qt , hence, 



sup 



3^% 



/ 1 



\Y(ty 



— sup 

t>o 



(1 



-at\k-l 



k>l 



< sup 

0<u<l 



k>l 



k'i 



,k-l 



k'i 



20 



L. LESKELA, PH. ROBERT, AND F. SIMATOS 



For < u < 1, the relation 



,k-l 



k>\ 



holds, hence 



1 -it 
— log u 

1 



o (1+*)" 

1 — 7} r oo 



dx, 



sup 

t>o 



Y(ty 



(x-logw)' 



< +oo. 



■ dx, 



The process (e~ at Y(t)) being a martingale, by convexity the process (e vat Y(t)~ n ) 
is a non-negative sub-martingale. For any r\ £ (0, 1) Doob's L p inequality gives the 
existence of a finite g(ry) > such that 

1 



Ei 



< q(v) sup 

t>0 



sup (e r " xt Y(t)- r ' 

,i>0 ^ 

The following result has therefore been proved 
Lemma A.l. For any < r\ < 1, 



;" Qt Ei 



< +00. 



Ei 



sup (>* (B^t) + I)"") 

t>0 v ' 



T = +00 



< +OO. 



On the event {r < +00}, hits a geometric number of times and then 

couples with a (p, A)-branching process conditioned on survival. On this event, 



sup 

t> 



p(e" Qt 0B CT (i) + l)-") 
v ' 

(e" Qt (fl ff (f) + I)"") ,sup (e" Qt (fl ff (f) + lp' 

<e" aT f 1 + sup (e" Qt + 1 



max [ sup 

^0<*<i 



where B' a (t) for i > r is the number of births in (r, t] of a (p, A)-branching process 
conditioned on survival and independent of the variable r, consequently 



E x 



sup 

t>0 



(e" at + I)"" 



r < +00 



<Ei (e" Q!T |T< +00) 



x i + : 



sup 

t>0 



(e" at (B a (t) + 1)-") 



T = + 0O 



In view of Lemma [A. 11 the proof of Proposition ^. 31 will be finished if one can prove 
that 

Ei (e" QT |r < +00) < +00, 
which actually comes from the following decomposition: under Pi( ■ \ t < +00), the 
random variable r can be written as 

l+G 

fc=l 

where G is a geometric random variable with parameter q = v j fj,z, (Tfe) is an i.i.d. 
sequence with the same distribution as the extinction time of a (p, A)-branching 
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process starting with one particle and conditioned on extinction and (E fJ/Zt k) are 
i.i.d. exponential random variables with parameter [iz- 

Since q is the probability of extinction of a {p, A)-branching process started with 
one particle, G + 1 represents the number of times (Z(t)) hits before going to 
infinity. This representation entails 

Ei (e" ar | t < +oo) = E (7(?7) G+1 ) where 7(77) = E (V'^+^W.i)) _ 

A (p, A)-branching process conditioned on extinction is actually a (1 — p, A)- 
branching process. See again Athreya and Ney [3]. Thus T\ satisfies the following 
recursive distributional equation: 

E x + l { t= 2 }{TiVT 2 ), 

where P(£ — 2) = 1 — p and E\ is an exponential random variable with parameter A. 
This equation yields 



P(Ti > t) < e~ At + 2A(1 - p) / P(Ti > t - u)e~ Au du 



[ P(Ti > i - 
Jo 



and Gronwall's Lemma applied to the function t (-> exp(Ai)P(Ti > t) gives that 

P(?i > t) < e (A_2A:p) * = e ( "~^ z)t 

hence for any < r\ < 1, 

Eite"" 7 ^ ' 



1-7? 

Since G is a geometric random variable with parameter q, E (j(r)) G ) is finite if and 
only if 7(77) < q. Since finally 

7fe) = ~ ^ (e" QTl ) < 7^ y 

one can easily check that 7(77) < q for 77 < rj*(v/^z) as defined by Equation ©, 
which concludes the proof of Proposition 



References 

[1] David Aldous and Jim Pitman, Tree-valued Markov chains derived from Galton- Watson 
processes, Annales de l'lnstitut Henri Poincare. Probabilites et Statistiqucs 34 (1998), no. 5, 
637-686. 

[2] Ceroid Alsmeyer, On the Galton- Watson Predator-Prey Process, Annals of Applied Proba- 
bility 3 (1993), no. 1, 198-211. 
[3] K. B. Athreya and P. E. Ney, Branching processes, Springer, 1972. 

[4] Maury Bramson, Stability of queueing networks, Lecture Notes in Mathematics, vol. 1950, 
Springer, Berlin, 2008, Lectures from the 36th Probability Summer School held in Saint-Flour, 
July 2-15, 2006. 

[5] Hong Chen and David D. Yao, Fundamentals of queueing networks, Springer- Verlag, New 
York, 2001, Performance, asymptotics, and optimization, Stochastic Modelling and Applied 
Probability. 

[6] T. D. Dang, R. Pereczes, and S. Molnar, Modeling the population of file- sharing peer-to-peer 
networks with branching processes, IEEE Symposium on Computers and Communications 
(ISCC'07) (Aveiro, Portugal), July 2007. 

[7] J. F. C. Kingman, The first birth problem for an age-dependent branching process., Annals 
of Probability 3 (1975), no. 5, 790-801. 

[8] John Lamperti, Continuous-state branching processes, Bulletin of the American Mathemati- 
cal Society 73 (1967), 382-386. 

[9] L. Lcskcla, Stochastic relations of random variables and processes, J. Theor. Probab. (2010), 
To appear. 



22 L. LESKELA, PH. ROBERT. AND F. SIMATOS 



[10 

[11 

[12; 

[13 
[14 

[15; 
[16 
[17; 
[18 

[19 



[20 
[21 



Laurent Massoulie and Milan Vojnovic, Coupon replication systems, Proceedings of SIG- 
METRICS'05 (Banff, Alberta, Canada), no. 1, June 2005, pp. 2-13. 

Olle Nerman, On the convergence of supercritical general ( C-M-J) branching processes, Z. 
Wahrschcinlichkeitstheorie verw. Gebiete 57 (1981), 365-395. 

J. Neveu, Erasing a branching tree, Advances in Applied Probability (1986), no. suppl., 101- 
108. 

R. Nunez-Qucija and B. J. Prabhu, Scaling laws for file dissemination in P2P networks with 
random contacts, Proceedings of fWQoS, 2008. 

Dongyu Qiu and R. Srikant, Modeling and performance analysis of bittorrent-like peer-to-peer 
networks, SIGCOMM '04: Proceedings of the 2004 conference on Applications, technologies, 
architectures, and protocols for computer communications (New York, NY, USA), ACM, 
2004, pp. 367-378. 

Philippe Robert, Stochastic networks and queues, Stochastic Modelling and Applied Proba- 
bility Scries, vol. 52, Springer, New- York, June 2003. 

L. C. G. Rogers and David Williams, Diffusions, Markov processes, and martingales. Vol. 2: 
ltd calculus, John Wiley & Sons Inc., New York, 1987. 

E. Seneta, On the supercritical branching process with immigration, Mathematical Biosciences 
7 (1970), 9-14. 

Florian Simatos, Philippe Robert, and Fabrice Guillcmin, A queueing system for modeling a 
file sharing principle, Proceedings of SIGMETRICS'08 (New York, NY, USA), ACM, 2008, 
pp. 181-192. 

Riikka Susitaival, Samuli Aalto, and Jorma Virtamo, Analyzing the dynamics and resource 
usage of P2P file sharing by a spatio-temporal model, International Workshop on P2P for 
High Performance Comptutational Sciences, 2006. 

David Williams, Probability with martingales, Cambridge University Press, 1991. 
Xiangying Yang and Gustavo de Veciana, Service capacity of peer to peer networks, Proceed- 
ings of IEEE Infocom'04, ACM, 2004, pp. 2242-2252. 



(L. Leskela) Aalto University, Department of Mathematics and Systems Analysis, PO 
Box 11100, 00076 Aalto, Finland 

E-mail address: lasse.leskelaOiki.fi 
URL: http://www.iki.fi/lsl 



(Ph. Robert, F. Simatos) INRIA Paris — Rocquencourt, Domaine de Voluceau, BP 105, 
78153 Le Chesnay, France. 

E-mail address: Philippe.RobertOinria.fr 
E-mail address: Florian.SimatosOinria.fr 
URL: http://www-rocq. inria.fr/~robert 



